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Abstract 

We study the large- N limit of a class of matrix models for dually weighted tri- 
angulated random surfaces using character expansion techniques. We show that for 
various choices of the weights of vertices of the dynamical triangulation the model can 
be solved by resumming the Itzykson-Di Francesco formula over congruence classes of 
Young tableau weights modulo three. From this we show that the large-iV limit im- 
plies a non-trivial correspondence with models of random surfaces weighted with only 
even coordination number vertices. We examine the critical behaviour and evaluation 
of observables and discuss their interrelationships in all models. We obtain explicit 
solutions of the model for simple choices of vertex weightings and use them to show 
how the matrix model reproduces features of the random surface sum. We also discuss 
some general properties of the large- iV character expansion approach as well as potential 
physical applications of our results. 
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1 Introduction 

The statistical mechanics of random surfaces has been of much interest over the years in 
the random geometry approach to two-dimensional quantum gravity and lower dimensional 
string theory. These models can be solved non-perturbatively by viewing discretized Riemann 
surfaces as Feynman graphs of N x N matrix models (see |l] for a review). The large- N limit 
of the matrix model exhibits phase transitions which correspond to the continuum limit of the 
dynamically triangulated random surface model and whose large- N expansion coincides with 
the genus expansion of the string theory. In this paper we will study a class of matrix models 
originally proposed by Das et al. |J in which the weighting of the coordination numbers 
of vertices (describing intrinsic curvature) of the dynamical triangulation of the surface can 
be varied. These "curvature" matrix models have been solved exactly in the large-iV limit 
for a Penner interaction potential by Chekhov and Makeenko ||, and more recently for an 
arbitrary potential by Kazakov, Staudacher and Wynter j|, [5] using group character expansion 
methods. 

There are many problems of interest which could be solved from the large- N limit of such 
matrix models. Foremost among these are questions related to the quantum geometry and 
fractal structure of the surfaces which contribute to the string partition function. For instance, 
the string theory ceases to make sense for string embedding dimensions D > 1. However, there 
is no obvious pathological reason for the existence of a barrier in the statistical mechanical 
model at D — 1, and it has been suggested || that the statistical mechanical model evolves 
into a different geometrical phase for D > 1, for instance a branched polymer (tree-like) 
phase, rather than the stringy (intrinsically two-dimensional) phase. It was argued in that 
curvature matrix models can represent such a transition by variation of the weighting of the 
vertex coordination numbers. The trees are to be thought of as connecting two-dimensional 
baby universes together in each of which the usual D < 1 behaviour is exhibited (see f?J for 
a review). Another problem is the exact solution of two-dimensional quantum gravity with 
higher curvature counterterms added to the action and the associated problem of the existence 
of a phase transition from the gravitational phase to a flat phase of the same string theory. 
By varying the vertex coordination numbers of a random lattice to a flat one, the genus zero 
problem can be obtained as the non-perturbative large- N solution of an appropriate curvature 
matrix model and it was demonstrated in || that there is no such phase transition. Besides 
these problems, there are many other interesting physical situations which require control 
over the types of surfaces that contribute to the random sum. 

The main subject of this paper is the curvature matrix model defined by the partition 
function 




(1.1) 
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where 

V 3 (XA) = \{XAf (1.2) 

and the integration is over the space of N x N Hermitian matrices X with A an invertible 
N x N external matrix. The Feynman diagram expansion of ( |1.1| ) can be written symbolically 

as 

^[A,*g=E n (v /2 t;f K) (i.3) 

where the sum is over all fat-graphs G3 made up of 3-point vertices, the weight associated to 
a vertex v* of the lattice dual to G 3 of coordination number q is 

t', = \^ , (") 

and M{Vq) is the number of such vertices in G3. This matrix model therefore assigns a weight 
t* whenever q 3-point vertices bound a face of the associated surface discretization, and it 
thus allows one to control the local intrinsic curvature R q = 7r(6 — q)/q of the dynamical 
triangulation. We will examine the structure of the solution at large-iV when these dual 
weights are arranged so that the only triangulations which contribute to the sum in ( |1.3|) are 
those whose vertices have coordination numbers which are multiples of three. 

There are several technical and physical reasons for studying such a model. The method 
of exact solution developed in [Q, [5] is based on a treatment of the large-iV limit of the 
Itzykson-Di Francesco character expansion formula |J . This formula is most naturally suited 
to matrix models of random surfaces whose vertices have even coordination numbers, and 
therefore the analysis in [||, 0, ^| was restricted to such situations. In the following we will 
show that the large-iV limit of the Itzykson-Di Francesco formula for the curvature matrix 
model fll-lfj can be used provided that one arranges the group character sum carefully. This 
arrangement reflects the discrete symmetries of the triangulation model that are not present 
in the models with vertices of only even coordination number. We shall see that this solution 
of the curvature matrix model actually implies a graph theoretical equivalence between the 
dynamically triangulated model and certain even coordination number models which were 
studied in [|], ||. We will also show how to map the Hermitian matrix model (|1.1|) onto a 
complex one whose character expansion is especially suited to deal with coordination numbers 
which are multiples of three and whose observables are in a one-to-one correspondence with 
those of the even coordination number matrix models. As a specific example, we solve the 
model explicitly for a simple power law variation of the vertex weights ( |1.4|) which agrees with 
expected results from the random surface sums and which provides explicit insights into the 
non-perturbative behaviour of observables of the curvature matrix model (|1 . 1|) , and also into 



the phase structure of the Itzykson-Zuber model ]I0 which is related to higher- dimensional 
Hermitian matrix models JTJ, [7], |TTJ. There are several physical situations which are best 
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described using such dynamical triangulations ||. The analysis which follows also applies to 
quite general odd coordination number models and establishes a non-trivial application of the 
techniques developed in @, 

The organization of this paper is as follows. In section 2 we derive the character expansion 
of the matrix model ( |1 . 1|) and discuss the structure of the Itzykson-Di Francesco formula in 
this case. We also show that certain symmetry assumptions that are made agree with the 
mapping of the model onto a complex curvature matrix model which demonstrates explicitly 
how the character sum should be taken. In section 3 we discuss the large-iV saddle point 
solutions of these matrix models. We show that the Itzykson-Di Francesco formula implies 
non-trivial correspondences between the Feynman graphs of the matrix model ( |1 . 1| ) and those 
of some models of random surfaces with even coordination numbers. We also establish this 
correspondence more directly using graph theory arguments. We analyse the critical behaviour 
of the matrix model when all vertices are weighted equally and discuss the various relations 
that exist among observables of the matrix models. We also demonstrate explicitly using the 
Wick expansion of ( |1 . 1| ) that the large- N solution of section 2 is indeed correct. In section 4 
we examine a simple situation where the dual vertices of the dynamical triangulation are not 
weighted equally. We show how this affects the large- iV saddle-point solution of the matrix 
model and thereby establish the validity of our solution for generic choices of the vertex 
coupling constants. Our results here may also be relevant to the study of phase transitions 
in the Itzykson-Zuber model ]7], |TTJ. In section 5 we briefly discuss the difficult problems 
which occur when dealing with matrix models that admit complex-valued saddle-points of 
the Itzykson-Di Francesco formula, and section 6 contains some concluding remarks and 
potential physical applications of our analysis. 



2 Large- N Character Expansion of the Partition Func- 
tion 

We begin by describing the character expansion method for solving the matrix model ( |1 . 1|) in 
the large- N limit j|, [5]. The external field A in (|1 . 1|) explicitly breaks the invariance of the 
model under unitary transformations X — > UXU> which diagonalize the Hermitian matrix 
X. Thus, unlike the more conventional matrix models with A = 1 for which the dual vertices 
v* are all weighted equally ]]]], the partition function ( |1 . 1| ) cannot be written as a statistical 
theory of the eigenvalues of X. The number of degrees of freedom can still, however, be 
reduced from iV 2 to N by expanding the invariant function e^ tT ( XA ) 3 i n characters of the 
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Lie group GL(N, C) as @] 

ef = c N £ X 3 [h) X{hi }{XA) (2.1) 

{hi} 

where denotes an irrelevant numerical constant and 



Ar/3 

n n (^ (£2) - ht ] 

0<ei<e 2 <2 ij=l 



(2.2) 



The sum in (|2.1| ) is over unitary irreducible representations of GL(N,C). They are charac- 
terized by their Young tableau weights {hi}f =1 which are increasing non-negative integers, 
< hi < h i+ i, and are defined by hi = i — 1 + b{ where 5, is the number of boxes in row % 
of the associated Young tableau. Because of the 3-valence coupling on the left-hand side of 
Q2.1J ), the sum is supported on those weights which can be factored into three groups of equal 



numbers of integers h\ e \ i = 1, . . . , y, where e = 0, 1, 2 labels their congruence classes modulo 
3. The GL(N, C) characters can be written using the Weyl character formula as 

X {hi) {Y) - A[y] (2.3) 
where jji G M. are the eigenvalues of the Hermitian matrix Y and 



A[y] = Y[( yi - Vj ) = detlyl- 1 ] (2.4) 

i<j 



is the Vandermonde determinant. 

Substituting (|2.1|) into and diagonalizing X, the integration over unitary degrees of 
freedom can be carried out explicitly using the Schur orthogonality relations for the characters. 
The resulting expression is a statistical mechanics model in Young tableau weight space which 
is a special case of the Itzykson-Di Francesco formula 0, || 

W N/ hh e - l)"/i°n / A \ -i^-^Z^iM+V) 

z H [\r q ] = c N £ ;t/ 2 ;, e } "\" x a [h) xm{A) A (2 . 5 ) 



The sum in ( J2.5[ ) is restricted to the even representations of GL(N, C), i.e. those with Young 
tableau weights which split into an equal number of even and odd integers h\ and h°, i = 
1, . . . , y. This restriction arises because of the Gaussian integration over the eigenvalues 
Xi G (— oo, oo) of the matrix X Q. Because of Q2.2|) , these weights must also split equally into 
their congruence classes modulo 3. The partition function (|2.5|) depends on only N degrees of 
freedom (the Young tableau weights hi) and thus the curvature matrix model (|1.1| ) is formally 
solvable in the large- iV limit. 

For definiteness, we consider the matrix model of dually weighted discretized surfaces for 
which the only graphs contributing in ( |L3p are those triangulations which have 3m, m = 
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1,2,..., nearest-neighbour sites to each vertex. This means that the vertex weights t* in QTT] 
are non-vanishing only when q is a multiple of three. To realize this explicitly in the matrix 
model ( |1.1| ), we take the external matrix A to be of the block form 

(A 1 / 3 \ 



A = A^ 



U3A 1 / 3 











■A 1 / 3 J 





V (J to 3 

where LU3 E Z 3 is a non-trivial cube root of unity and A 1 / 3 is an invertible A x/ N 
matrix. The character ( |2.3| ) can then be evaluated explicitly to get [Q 



(2.6) 



3 x y Hermitian 



and the statistical sum (E.5I) becomes 



(A) sgn 



AT/3 



n n (* 

0<ei<e 2 <2 i,j=l 



(^2) 



/A[ft«>] X'r^l W)\ 



e=0,l,2 



x e^i 



EiM|log(&)+£log(3)] 



(2.7) 



(2.8) 



Notice that the sign factors from fl2.2| ) and (277) cancel each other out in 

To treat the large- N limit, we assume that at iV = 00 the sum over representations in 
(|2.8|) becomes dominated by a single, most-probable Young tableau {hi}. The problem that 
immediately arises is that the groupings of the Young tableau weights in ( |2.8|) into equal 
numbers of even and odd integers and into equal numbers of mod 3 congruence elements 
need not occur symmetrically. There does not appear to be any canonical way to split the 
weights up and distribute them in the appropriate way. However, it is natural to assume, 
by the symmetry of the weight distributions in ( |2.8|) , that the saddle-point localizes around 
that configuration with equal numbers of even and odd weights, each set of which groups 
with equal numbers into their mod 3 congruence classes hi and h° , i — 1, . . . , y. It is 
also natural, by symmetry again, to further assume that although the different sets of weights 
h\ , h°^ do not factorize and decouple from each other, the groupings of the weights into 
even and odd mod 3 congruence classes are distributed in the same way. We therefore assume 
this symmetrical splitting in (|2.8| ) and write the statistical sum over the mod 3 congruence 
classes of weights h[ . We also make the additional assumption that the product of weights 
from two different congruence class groupings contributes the same in the large-iV limit as 
the Vandermonde determinant for a single grouping, i.e. that 



n(X ( 



h) 



(«*) 



hi 



n W ei) - $ 



(*i) 



(2.9) 
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for £2 7^ t\. As we shall see, these symmetrical grouping assumptions are indeed valid and 
lead to the appropriate solution of the curvature matrix model ( |1 . 1| ) at large- N. 

We now rescale hi — > iV • hi in (|2.8| ), simplify the factorials for large- N using the Stirling 
approximations h\\ ~ gM^s^- 1 )/ 2 anc i /^l ^ e (/i+ 5 )logft.-ft,^ an( j a pply the above symmetry 

assumption retaining only the leading planar (order 1/N) contributions to (|2.8| ). After some 
algebra, the partition function (|2.8| ) for the symmetrically distributed weights h^ in the 
large- N limit can be written as 

4°) „ Yl e"* 3 "^ (2.10) 

where the effective action is 

i q AT/3 -, JV/3 

^[/> (0) ] = --log A + £ log^f - + -L £ fcf) [2 log A + log 

i<jf j=l 



: , A 



3 , r 



and we have introduced the Itzykson-Zuber integral 

.La^w'-v a[q] 



(2.11) 



^'^U df7 ^ *™ = W^spTO (2 ' 12) 

where a, are the eigenvalues of the matrix A and a j = log Oj . In arriving at ( [2.10|) we have 
used the Vandermonde determinant decomposition 

AT/6 

A[/i (£) ] = A[/i e(£) ]A[/i o(e) ] ]J (hf t] - hf e) ) (2.13) 

and ignored irrelevant terms which are independent of A and the h's. 

In the case of the even coordination number models studied in [5|, H the natural weight 
distribution to sum over in the Itzykson-Di Francesco formula ( |2.5| ) is that of the original 
even representation restriction (see subsection 3.1). In our case the natural weight distribu- 
tion appears to be that of the mod 3 congruence classes. However, in contrast to the even 
coordination number models, there is no way to justify this from the onset. We shall see 
later on that the appropriate distribution for the Itzykson-Di Francesco formula at large- iV is 
essentially determined by the discrete symmetries of the given curvature matrix model. Some 
evidence for the validity of this assumption comes from comparing the Hermitian model here 
with the complex curvature matrix model 

Z C [XX] = (2n\)- N2 [ d0 dtf e~T tr H^Bf (214) 



where the integration is now over the space of N x N complex-valued matrices <f), and B 
is another external N x N invertible matrix. Again one can expand the invariant, cubic 
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interation term as in ( |2.1| ) and the character expansion of ( |2.14 ) is |T2"f 

.ijV(JV-i)+£./H 

N 



Zc[XX 



h={h(*)} 



Y\ N h I / A N 

a -' X{h}(A) X{h }(B) X 3 [h] [ - 



A[h] 



(2.15) 



Now the Gaussian integration over the eigenvalues of the positive definite Hermitian matrix 
(f)'4> goes over [0, oo) and so there is no restriction to even and odd weights in ( j2.15| ), only the 
restriction to mod 3 congruence classes because of the 3- valence coupling. If we take B = 1, 



so that X{h}(B) oc A[h], and use the decomposition 
becomes 



then the character expansion (|2.15 ) 



N 



z c [x,t;]=c N x-^ N ^ £ nwn- 

/j={hW}i=l 6=0,1,2 J] 



A[h^ e 



X 



(2.16) 



Thus in this case the congruence classes of weights completely factorize, and we can now 
naturally assume that the N = oo configuration of weights is distributed equally into these 
three classes. In this sense, the complex matrix model (|2.14j) is a better representative of the 
dually-weighted triangulated random surface sum, because one need not make any ad-hoc 
assumptions about the weight distribution. The character expansion Q2.16j ) suggests that the 
correct statistical distribution for the dynamical triangulation model at large- N is over hf\ 
as was assumed above, and we shall discuss the explicit relationship between these Hermitian 
and complex matrix models in the next section. We now rescale hi — > N ■ hi and simplify 
(|2.16| ) for large-iV as for the Hermitian model. After some algebra, the partition function 
(|2.16| ) for the symmetrically distributed weights h^ in the large- limit can be written as 



7(0) 

J c 



him 



N 2 Sc[h( 0) ] 



(2.17) 



where the effective action is 



Sc[h {0) ] 
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N/3 



N / 3 

2 , (o) 



Kj 1=1 



log (X^h 



(0) 



+ A^ l0g/ 



Nh {0) 



A 



(2.18) 



3 Saddle-point Solutions 



If the characteristic h's which contribute to the sum in fl2.10| ) are of order 1, then the large- A" 



limit of the Itzykson-Di Francesco formula can be found using the saddle-point approximation. 
The attraction or repulsion of the A-dependent terms in fl2.8| ) is compensated by the Vander- 
monde determinant factors, so that, in spite of the unsymmetrical grouping of Young tableau 
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weights that occurs, the partition function is particularly well-suited to a usual saddle-point 
analysis in the large- N limit. In this section we shall see that the symmetry assumptions made 
above give a well-posed solution of the matrix model at iV = oo. The partition functions ( |1 . 1| ) 
and ( |2.14j ) are dominated at large- N by their saddle points which are the extrema of the ac- 
tions in ( [2. 11| ) and fl2.18| ). In what follows we shall determine the saddle-point solutions of 
both the Hermitian and complex matrix models and use them to examine various properties 
of the random surface sum. 



3.1 Hermitian Model 

To find the saddle-point solution of the Hermitian matrix model, we minimize the action 



( 2.11 ) with respect to the Young tableau weights. Then the stationary condition J^jgj- = 



leads to the saddle-point equation 



o N/3 -, 

^T, W ^ m = -log(A 3 /4 0) ) -2^1 0) ) (3.1: 

i 3 



where we have introduced the Itzykson-Zuber correlator 



„ m = (3 . 2) 

oh) 



To solve ( |3.1|) , we assume that at large- N the Young tableau weights h 



(0) , TV 

' L i ■ ■ ■ i 3 j 



become distributed on a finite interval h G [0, a] on the real line and introduce the normalized 
spectral density pu(h) = ^^y, where h(x) G [0, a] is a non- decreasing differentiable function 



of x G [0, 1] with h(3i/N) = N ■ hf ] . Discrete sums over h) are then replaced with integrals 



over h by the rule -j| Ylj=i ~ > Jo dh. Notice that since the h^s are increasing integers, we have 
< Puih) < 1, and so the spectral distribution is trivial, pn(h) = 1, on some sub- interval 



[0,6] with < b < 1 < a [13|. The saddle-point equation ( |3.1|) then becomes an integral 



equation for the spectral density, 

£ dU = - log(A3/l) - log {h~b) - 2:FW ' he[b > a] (3 - 3) 

where we have saturated the spectral distribution function at its maximum value pu{h) = 1 
on [0, b\. The saddle-point solution of the matrix model can thus be determined as the solution 
of the Riemann-Hilbert problem for the usual resolvent function 
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which is analytic everywhere in the complex /i-plane away from the support interval [0, a] of 
Ph where it has a branch cut. The discontinuity of the resolvent across this cut determines 
the spectral density by 

, PH(h') 

'0 



H H (h±iO) = J* dti ^—± T mp H (h) , he[0,a] (3.5) 
jo h — h' ' 



In contrast to the more conventional Hermitian one-matrix models the Riemann- 
Hilbert problem (|3.3| ) involves the unknown Itzykson-Zuber correlator T{h) which must be 
determined separately in the large- N limit. As shown in j|, 0], it is determined by the vertex 
couplings (Of) through the contour integral 



9? = ^ = = J £ ^ e q(HH{h)+Hh)) , Q>1 (3-6) 

where the closed contour C encircles the support of the spectral function pn{h), i.e. the cut 
singularity of 7i#(/i), with counterclockwise orientation in the complex /i-plane. Note that 
( |3.6j ) is the large- N limit of a set of weights of size N/3 and it follows from the identity Q 

tr A* = E "J fjy- (3-7) 

fc =i X{M0)/3}1 A J 

where 

(4 0) (?))i = ^ 0) + 3g^ (3.8) 

In the next section we shall discuss the evaluation of T(h) and the corresponding structure of 
the curvature matrix model as the vertex weights t* are varied. Notice that, strictly speaking, 
in most cases of interest the characters corresponding to a specific choice of vertex weightings 
cannot be represented via matrix traces as in (1.4) and need to be defined by an analytical 



continuation. This can be accomplished by using the Schur-Weyl duality theorem to represent 
the GL(N, C) characters as the generalized Schur functions 

X{h}[t*)=det[P hk+1 „ i [t*)) (3.9) 
where P n [t*] are the Schur polynomials defined by 

(oo \ oo 

Nj2z%) =Y, znp nlt*} (3-10) 
g=l / n=0 

When the weights t* are given by (|1.4j ), the Schur functions ( |3.9| ) coincide with the Weyl 



characters (2.3). 



Our first observation here is that the saddle-point equation (|3.3j) is identical to that of 



the even-even coordination number model discussed in [0L K| which is defined by replacing the 
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cubic interaction matrix potential V${XA) in by 



V even (XA) = J2t2 q (XA^ 

9=1 



where 



A 1 ' 2 o 
o -A 1 ' 2 



(3.11) 



(3.12) 



with A 1 ! 2 an invertible y x Y Hermitian matrix, and 



tr ~ 



2qt 2q = 2qt* 2q = 2-A« 



N 



(3.13) 



The curvature matrix model with potential ( |3.11| ) generates a sum over random surfaces where 
arbitrary even coordination number vertices of the primary and dual lattices are permitted 
and are weighted identically, so that the model Q3.11|) is self-dual. In fact, the same saddle- 
point equation arises from the interaction potential 



V A (XA) = ^(XA {2) Y 



(3.14) 



The matrix model with potential V4 generates discretizations built up from even-sided poly- 
gons with 4-point vertices. 

As shown in J3J, the Itzykson-Di Francesco formula for the curvature matrix model with 
potential Q3.11| ) follows from replacing X 3 [h] in ( |2 . 5|) by 



X{h}(AV>) = X {^j(A) X{^i}(i) sgn 



' N/2 

n 0* 



(3.15) 



Now the even and odd weights can be naturally assumed to distribute equally, and the effective 
action is 



*Seven [h 



, 2 N/2 1 N/2 

■i^g\+—Y: \o g (ht hi [iog(xht) - 1] +4 log j 



Nh e 



A 



(3.16) 



Defining a distribution function for the N/2 weights hf analogously to that above and varying 
the action ( |3.16[) yields the saddle-point equation (|3.3f) with A 3 — > A. In this case the Itzykson- 
Zuber correlator J-{h) is determined as in ( |3.6|) but now with qt* defined to be equal to (|3.13|) 
(i.e. the large- N limit of a set of weights of size N/2). 

Likewise, in the case of the 4-point model ( |3.14j ), we replace X s [h] by 



X A [h) 



sgn 



AT/4 

n n (^ 2) - 

0<Mi<M2<3 i,j=l 



(3.17) 
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where now the character sum is supported on weights that factor equally into their congruence 
classes h!f\ p, — 0, . . . , 4, i — 1, . . . , modulo 4. Because of the original weight constraint 
of the Itzykson-Di Francesco formula, this means that the even weights hf distribute equally 
into the mod 4 congruence classes p = 0, 2 and the odd weights h° into the classes p, — 1, 3. 
Again, assuming these weights all distribute equally leads to the effective action 



i i N/2 1 N/2 

S,[h e ] = --l gX+—J2^g(ht-^) + —J2ht [MA 2 ht) -l]+21ogJ 

i<j i=l 



Nh e ~ 
— , A 



(3.18) 



Introducing a distribution function for the N/2 weights hi again leads to precisely the same 
saddle-point equation ( |3.3p with A 3 — > A 2 . Here and in the even-even model above, unlike the 
3-point model of the previous section, the large- N configuration of weights naturally localizes 
onto h\. 

These three matrix models therefore all possess the same solution at iV = oo, i.e. their 
random surface ensembles of genus zero graphs are identical. Their 1/N corrections will differ 
somewhat because of the different ways that the weights split into the respective congruence 
classes in the three separate cases. The genus zero free energies are related in a simple 
fashion. From (|2.11|) and the definition of the Hermitian distribution function, the genus zero 
free energy for the 3-point model is 

S H [X, £) = l -( l - ( a Tdh dti PH {h) PH {h') log \h - h'\ + f dh p H (h)h [log(A 3 /i) - 1 
o \2 Jo J o Jo L 

+2j Q a dh p H (h) log I c [h, A]) - ^ log A 

(3.19) 

where pn{h) is the solution of the saddle-point equation ( |3.3j ) and log I c [h, A] = 3 2 \ogI[ Nh ^ 0) , A] 
is the Itzykson-Zuber integral (|2.12|) at = oo. Similarly, the genus zero free energies for 
the even-even and 4-point models can be written using the same distribution function pu{h). 
Now, however, the rules for replacing sums by integrals differ. For the even-even and 4- 
point models the normalized sums corresponding to a spectral density normalized to unity 
are jjY^ili an d the N = oo Itzykson-Zuber integral is log I c [h, A] = 2 2 log I[^f-, A]. Taking 
this into account along with the change of A in the three cases, we find that the free energies 
of the three models discussed above are all related by 

3SW[A 2 / 3 , t* q ) = 2S 4 [\, t* q ] = S cvcn [X 2 , t* q ] (3.20) 

It is quite intriguing that the Itzykson-Di Francesco formula naturally implies these re- 
lations by dictating the manner in which the Young tableau weights should decompose with 
respect to the appropriate congruence classes in each case. This is reflected in both the overall 
numerical coefficients and the overall powers of A that appear in ( |3.20| ). In the next subsection 
we shall give purely graph-theoretical proofs of these relationships. This yields a non-trivial 
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Figure 1: A 3-colouring of a graph in M3 (thick lines) and its associated dual graph in M$ (thin 
lines). 



verification of the assumptions made in section 2 about the splitting of the Young tableau 
weights for the dynamical triangulation model. 



3.2 Graph Theoretical Relationships 

Before pursuing further the properties of the curvature matrix model above, we present a 
direct, graphical proof of the relationship (|3.20| ) between the generating functions for the 



3-point, 4-point, and even-even planar graphs. We denote the ensembles of spherical graphs 
in these three cases by A4 3 , M.^ and A4 even , respectively. We first discuss the mapping 
— > M.4. Consider a planar graph G G Ai%. Using a version of the colouring theorem for a 
spherical topology, we can colour the lines of G by three colours, R, B, and G, to give a planar 
graph G c with labelled lines (Fig. 1). Because G consists of polygons whose numbers of sides 
are all multiples of three, the colouring can be chosen so that the three colours always occur 
in the sequence (R,B,G) as one goes around any polygon of G c with clockwise orientation. We 
can now contract the two 3-point vertices bounding each R-link of G c into a 4-point vertex 
(Fig. 2). A polygon of 3m sides with 3-point vertices in G c then becomes a polygon of 2m 
sides with 4-point vertices. The resulting 2-coloured graph G c thus belongs to M.4. Notice 
that if each link has a weight A associated to it, then the effect of the R-contractions is to 
map A 3 — > A 2 on M3 — > M4. 

Conversely, suppose G e .M4. We can form a 2-coloured graph G c , with colours B and G, 
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Figure 2: A 2-colouring of a graph in M.4 (thick lines), its associated dual graph in M% (thin lines), 
and the corresponding diagram in -M e ven (dashed lines). The 4-point vertices denoted by solid 
circles are contracted from the R-links of the corresponding graph of M.3 depicted in Fig. 1. 



such that the colours alternate along the lines of G c (Fig. 2). There are two possible ways 
to orient lines of G c , by denning a direction to them from either B to G or from G to B at 
each vertex. The 4-point vertices of the resulting oriented graph can then be split into a line, 
labelled by R, bounded by two 3-point vertices. There are two ways of doing this, by splitting 
the 4-point vertex either vertically or horizontally with respect to the given orientation (Fig. 
3). Thus to each 2-coloured graph G c G M.^ there corresponds two distinct (topologically 
inequivalent) 3-coloured graphs G c G M.3. There are also three distinct 2-coloured graphs 
G c G M.4 for each 3-coloured graph G c G M 3 corresponding to the three possible choices of 
contraction colour R, B or G. This therefore defines a three-to-two mapping on .M4 — > 
and is just the statement of the first equality of ( |3.20|) . 

Actually, there exists a simpler, line mapping between these two ensembles of graphs in 
terms of dual lattices. Let .M3 denote the collection of planar graphs dual to those of M.3 (i.e. 
lattices built up from triangles that form 3m- valence vertices), and M\ the dual ensemble 
to A^4 (i.e. the graphs formed of squares that meet to form vertices of even coordination 
number). The sets M.\ and M.\ are generated, respectively, by the curvature matrix models 
with matrix potentials 

00 00 
vt\XA) = Y,tl q {XA^ , vt\XA)=J2t2 q (XA 4 )^ (3.21) 

q=l g=l 

where the N xN matrix A m is denned by jjA*^ = 5^. The corresponding generating functions 
S$\\,i*] and Si* ] [\,t* q } coincide, respectively, with (|XL| ) and (pT8|) [§. A gain the lines of 
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R ^ — 





B 



Figure 3: The two possible 4-point vertex splittings which respect the given B and G colour orien- 
tation. 



the graphs of Ai% can be 3-coloured, and the deletion of all R-coloured lines gives a map 
Ai^ —>■ Ail (Figs. 1 and 2). This mapping is three-to-one because of the three possible 
choices of deletion colour. The inverse map Ail — > Ai\, defined by inserting an R-link across 
the diagonal of each square of the graphs of Ail, is then two-to-one because of the two possible 
choices of diagonal of a square. 

The correspondence between the ensembles of graphs Ai^ and A4 even is similar and has 
been noted in || (see Fig. 2). Given a graph G* G Ail, we choose a diagonal of each square 
of G*. Connecting these diagonals together produces a graph in -M e ven- Equivalently, we can 
place vertices in the face centers of the corresponding dual graph G G Ai^ and connect them 
together by drawing lines through each of the 4-point vertices of G (so that A 4 — > A 8 from the 
splitting of these 4-point vertices). In this way we obtain a map Ai^, Ail — > Ai evcn where the 
vertices and face centers of graphs of -M e ven correspond to the vertices of graphs of Ail, or 
equivalently the faces of Ai^. Because there are two distinct ways of choosing the diagonal 
of a square of G* G Ail (equivalently two ways of splitting a 4-point vertex of G G Ai± 
with respect to a 2-colour orientation analogously to that in Fig. 3), this mapping is two-to- 
one and corresponds to the second equality of ( p.20[ ). In particular, there is a three-to-one 
correspondence between graphs of A^A^ and .Meven- 



We stress that these graphical mappings are only valid on the sphere. In terms of the 
Itzykson-Di Francesco formula, this means that the 0(1/N) corrections to the large- N saddle- 
point solutions of the corresponding curvature matrix models will differ. These non-trivial 
correspondences are predicted from the matrix model formulations, because in all cases we 
obtain the same N = 00 saddle-point equations but different splittings of Young tableau 
weights into congruence classes thus leading to different overall combinatorical factors in front 
of the graph generating function. Note that in the case of the 4-point and even-even models 
the vertex weights t* are mapped into each other under the above graphical correspondence 
because of the simple line map that exists between -M e ven and Ail- In the case of the mappings 
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onto M. 3 the weights t* defined in ( |3.6| ) are mapped onto those defined by ( |3.13| ) because of 
the contraction of 3m-sided polygons into 2m-sided polygons. 



3.3 Complex Model 

As mentioned in section 2, it is useful to compare the Hermitian 3-point curvature matrix 
model with the complex one since in the latter case the Young tableau weights completely 
factorize and the splitting of them into mod 3 congruence classes appears symmetrically. The 
saddle-point solution in the case of the complex curvature matrix model is identical in most 
respects to that of the Hermitian models above. Now the spectral density pc{h) for the Young 
tableau weights hf \ i = 1, . . . , ^, obeys the saddle-point equation 



i 



b dh ' = - MA 3/2 />) - log (^-^J - \Hh) , he[b, a] (3.22) 

where the logarithmic derivative J-{h) of the Itzykson-Zuber integral is defined just as in 
(|3.2j ). The solution for J-'(h) = is thus identical to those above. In particular, working out 
the free energy as before we see that 

S C [X, t* q ] = AS H [p c ; Vx, Q - £ dh p c (h) log I c [h, A] (3.23) 



The combinatorical factors appearing in (|3.23|) can be understood from the Wick ex- 
pansions of the Hermitian and complex curvature matrix models. First consider the case 



JF(/i) = 0. One factor of 2 appears in front of Sh in ( 3.23 ) because the number of inde- 
pendent degrees of freedom of the N x N complex matrix model is twice that of the N x N 
Hermitian curvature matrix model. The other factor of 2 arises from the mapping of Feynman 
graphs of the complex matrix model onto those of the Hermitian model (Fig. 4). At each 
6-point vertex of a complex graph we can place a $ line beside a line to give a graph with 
3-point vertices and "thick" lines (each thick line associated with a </>'</> pair of lines). This 
maps the propagator weights as A 6 — > A 3 and there are 3!/3 = 2 distinct ways of producing 
a complex graph by thickening and splitting the lines of a graph in Ai 3 in this way. This is 
the relation ( [3.23D for J c = 1. In the general case, there is a relative factor of 4 in front of 



the Itzykson-Zuber correlator in ( |3.22|) because the Feynman rules for the complex curvature 
matrix model associate the weights 

«? = ;f KT p- 24 > 

to the dual vertices v* of the Feynman graphs (compare with ( |1.4j )). Thus for qt* 7^ 1, the 
free energy of the complex model differs from ( |3.19| ) in a factor of 4 in front of the integral 
involving the Itzykson-Zuber integral I c [h, A]. 
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Figure 4: The mapping of a complex model 6-point vertex onto a Hermitian model 3-point vertex. 
The incoming lines represent the <jy fields, the outgoing lines the <f> fields, and the solid circle 
denotes an insertion of the matrix A 1 / 2 . 

The complex curvature matrix model with the above rules thus generates "checkered" 
Riemann surfaces corresponding to the double-lined triangulations shown in Fig. 4. A 3- 
colouring of a graph of is now described by a sort of 6-co louring whereby each colour 
R, B, and G is assigned an incoming and outgoing orientation at each vertex. Again these 
relations are only valid at A = oo. However, the models lie in the same universality class (as 
with the even coordination number models discussed above) and consequently their double 
scaling limits will be the same. In particular, the equivalences (|3.20|) and (|3.23|) will hold to 
all genera near the critical point fl4| (see the next subsection). 

3.4 Critical Behaviour and Correlation Functions 

We will now discuss to what extent the above models are universal. First, let us find the 
explicit solution of the saddle-point equation (|3.3|) in the case J-'(h) = 0] (i.e. A — A — 1 
and qt* = 1). Notice that in this case the dual potentials in (|3.11| ) and ( |3.21| ) become the 
Penner-type potentials J3], |9| 

V even (XA) = -~ log (l - {XA^f) , V£\XA) = -\ log (l - (XA 4 ) 2 ) 

1 1 (3.25) 

vt\XA) = --log(l-(XA 3 f) 

so that the above arguments imply the equivalences between the curvature matrix model for 
a dynamical triangulation and those of Penner models. The behaviour of the saddle-point 
solution when the dual vertices v* are not weighted equally will be examined in the next 
section. 

The saddle-point equation ( |3.3D determines the continuous part of the resolvent function 
Tinih) across its cut. Using the normalization J* a dh pn{h) = 1 of the spectral density in 
(|3.4j ) implies the asymptotic boundary condition Tinih) ~ l//t at |/t| — > oo. It follows that 
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the solution for the resolvent function is given by |1| 



n„(h) 



log 
log 



h 



h-b / 

' h x 
h-b 



ds 



(h-a)(h-b) 



+ 



c 2-ai s — h\ (s — a)(s — b) Jb 
1 



-Tdh 
Jb 



, PH{h') 

-h! 



ds 



(3.26) 



c 2ni s — h\ (s — a)(s — b) 



(h - a)(h - b) f , . o . . / s 
J log(A 3 s) + log 



s-b 



The contour integrations in (|3.26|) can be evaluated by blowing up the contour C to infinity 
and catching the contributions from the cuts of the two logarithms, on (— oo, 0] and [0,b], 
respectively, and also from the simple pole at s — h. Keeping careful track of the signs of the 
square roots on each cut and taking discontinuities across these cuts, after some algebra we 
arrive at 



n H {h) 



log(A 3 /i) + y/(h-a)(h-b) 



dx 



log 



(h+ Vab + ^J(h - a)(h - 6))" 



{a — b 
A 3 (v^ + Vbf h (h(a + b) - 2ab + 2y/ab(h - a)(h - b) 



(3.27) 



It remains to solve for the endpoints a, b of the support of the spectral distribution function 
in terms of the parameters of the matrix potential. They are determined by expanding the 
function (|3.27 ) for large-/;., and requiring that the constant term vanish and that the coefficient 
of the 1/h term be 1. This leads to a pair of equations for a, b 



3AV = t]-1 



where we have introduced the positive endpoint parameters 

1/^/^2 1 



Vb 



v 



4 



(3.28) 



(3.29) 



The cubic equation for r] in ( 3.28 ) can be solved explicitly. The original matrix model free 
energy is analytic around A = 0. We shall see below that it is analytic in rj, so that we should 
choose the branch of the cubic equation in ( |3.28| ) which is regular at A = 0. This solution is 
found to be 

iVz (V/ 3 



V 



. V — 



gA 6 /^ 



(3.30) 



where 



' + JtttV / 81A 6 -4 



(3.31) 



6A 6 54A 9 

The solution ( |3.30| ) leads to the appropriate boundary condition r/(A = 0) = 1, £(A = 0) = 0, 
i.e. a = b = 1, as expected since in this case the spectral density should be trivial, puih) = 1. 
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The solution ( |3.30| ) is real- valued for A < A c = (2/9) 1 / 3 which is the regime wherein all 
three roots of the cubic equation for rj in ( |3.28| ) are real. For A > A c there is a unique real- 
valued solution which does not connect continuously with the real solution that is regular 
at A = 0. At A = A c the pair of complex-conjugate roots for A > A c coalesce and become 
real-valued. Thus for A > A c the one-cut solution above for TCh is no longer valid and the 
spectral density pn becomes supported on more than one interval. The value A = A c is 
therefore identified as the critical point of a phase transition of the random surface theory 
which corresponds to a change of analytic structure of the saddle-point solution of the matrix 
model. As usual, at this critical point the continuum limit of the discretized random surface 
model, describing the physical string theory, is reached. To determine the precise nature of 
the geometry that is obtained in this continuum limit, we need to study the scaling behaviour 
of the matrix model free energy near the critical point. In the following we shall present two 
calculations illustrating this critical behaviour. 



3.4.1 Propagator 

As shown in ||, |5j, in the general case the resolvent function TCii(h) and the Itzykson-Zuber 
correlator ^(h) completely determine the solution of the random matrix model in the large- 
N limit. For instance, expectation values of the operators jjX 2q can be obtained from the 
contour integral 

/ ^ X 2q ) =^-t^Lh« e«"*W (3.32) 
\N / Ma q Jc 2ni v ' 

where the normalized average is with respect to the partition function (|1.1|) and the decompo- 
sition ( p.6| ). The expression ( |3.32j ) holds only for the even moments of the matrix X because it 
is only for such moments that the character expansion of the right-hand side of ( |3.32p admits 



a symmetrical even-odd decomposition as in the Itzykson-Di Francesco formula M. Although 
a character expansion for the odd moments can be similarly derived, the weights will not 
split symmetrically and there does not appear to be a general expression for them in terms 
of the weight distribution function at large-A^. For the even coordination number models 
above, the reflection symmetry X — > —X of the matrix model forces the vanishing of all odd 



moments automatically and (|3.32| ) represents the general expression for the complete set of 
observables which do not involve the external field A at large- N . In terms of the character 
expansion formula, the character sum for (^X 2q+1 ) vanishes in the case of an even potential 
because of the localization onto the even representations of GL(N, C), whereas in the 3-point 
model (jjX 2q+1 ) m 3 7^ because its character sum is supported on mod 3 congruence class 
representations. 

The contour integral representation ( |3.32| ), with Tin replaced by Tic also holds for the 



correlators (jj((p^(p) q )c of the complex curvature matrix model. These averages form the 



18 



complete set of A-independent observables at large-iV in this case because of the charge 
conjugation symmetry <fi — > $ of the complex model (which is the analogue here of the 
reflection symmetry of the even coordination number models). Again this indicates that the 
complex curvature matrix model is better suited to represent the triangulated random surface 
sum, because it is more amenable to an explicit solution in the large- N limit. Furthermore, its 
complete set of correlation functions coincide at N = oo with those of the equivalent matrix 
models defined above with even potentials, after redefining the vertex weights according to 
(0) and (pi). 

The contour integral on the right-hand side of ( |3.32| ) can be evaluated by blowing up the 
contour C, expanding the resolvent function (ft.27|) for large-/;., and computing the residue at 
h = oo. For example, for q = 1 we have 

C4X 2 ) =x(U(h)) (3.33) 



\N /m 3 V2 

where the weight average (h) = (jfj^i^ih^) = Jq dh pn{h)h corresponds to the coefficient 
of the 1/h 2 term in the asymptotic expansion of ( |3.27| ). This result also follows directly from 
differentiating the free energies 

= 2*1^*3 , (£*V) o = A"!*,[A,*g (3-34) 



and using the saddle-point equations. Notice that ( |3.34| ) and the free energy relations ( |3.20| ) 
and ( |3.23| ) imply, in particular, that the matrix propagators of the 3-point, 4-point, even- 
even and complex matrix models coincide with the appropriate redefinitions of the coupling 
constants A and t*. However, the equivalences of generic observables in the various models 
will not necessarily hold. We shall return to this point shortly. 

Using ( |3.271) and (|3.28p , the propagator (|3.33p can be written as 



\N /Ms V 3 3 

Expanding (|3.35|) using ( |3.30|) as a power series in A yields 

/ ^X 2 ) = A + A 7 + 6A 13 + 54A 19 + 594A 25 + £>(A 31 ) (3.36) 

\N I M-i 

We have verified, by an explicit Wick expansion of the matrix propagator (up to and including 
the order shown in ( |3.36| )) using the 3m-sided polygon constraints on the matrix traces fll.4|) , 
that ( p.36| ) indeed coincides with perturbative expansion of the curvature matrix model (|1.1| ) 
and thus correctly counts the planar 3-point fat graphs consisting of only 3m-sided polygons. 
It also agrees with the Wick expansions of the even coordination number models above M 
and of the complex curvature matrix model. 
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To examine the scaling behaviour of observables of the matrix model near the critical 
point, we introduce a renormalized, continuum cosmo logical constant A and momentum II by 



A 6 = A*(1-A) 



v = vc- n/2 



(3.37) 



where rj c = r){\ = A c ) = 3/2. We then approach the critical point along the line A, 11(A) — ► 0, 
where the function 11(A) is found by substituting the definitions fl3.37| ) into ( |3.3(J| ) to get 



n(A) 



5l/3+ 9(A. + l)_ 6 + ^/ 5l/3 _9(A. + l)- 



with 



S = 27 A* + 1 + J-Al - 2A 2 - A, 



A* 



A 



1 - A 



(3.38) 



(3.39) 



Substituting ( |3.37[) and (|3.38|) into (|3.35|) and expanding around A = yields after some 
algebra 

/ ^x 2 \ = -A 3/2 + . . . = c • A 1 "** + . . . (3.40) 
\N I Mi 9 

where the dots denote terms which are less singular as A — > 0. The leading non- analytic 
behaviour in ( 3.40|) identifies the critical string exponent of this random surface model as 



7str = "1/2 



(3.41) 



so that the system in the continuum limit represents pure two-dimensional quantum gravity 
p] . The same argument also applies to the 4-point and even-even models with the appropriate 



redefinitions of cosmological constants in ( |3.37|) . Thus the equal weighting of all vertices in 
these random surface models leads to curvature matrix models in the same universality class 
as the more conventional Hermitian one-matrix models of two-dimensional quantum gravity. 



This agrees with recent numerical results in [15 



3.4.2 General Correlators 

We now turn to a more general discussion of the evaluation of observables in the matrix 
models above. The relations between correlators involving the external matrix A are a bit 



more subtle than those represented by the contour integrations in ( |3.32|) which coincide in all 
four matrix models for all q. Because of the simple line maps that exist between the ensembles 
A^3, and .M e ven, all expectation values are the same in these models, i.e. 

Analytically, this equality follows from the fact that these correlators are given by derivatives 
of the free energies of the matrix models with respect to the weights t*. The powers in (|3.42|) 
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can be understood from the fact that the reflection symmetry — > —A^ of the .M4 and 
A4even models (which restrict the non-zero averages to (^(X A^) 2q ) M4jMcvcn ) corresponds 
to the Z 3 -symmetry A^ — > w 3j A®, CU3 G Z 3 , of the .M 3 model (which restricts the non- 
vanishing observables to (jj(XA^) 3q ) M3 and (^(0V^^) 39 )c)- In terms of the coloured 
graphs of subsection 3.2, these discrete symmetries correspond to permutations of the colours 
of graphs in each ensemble. It is these symmetries that are ultimately responsible for the 
localization of the Itzykson-Di Francesco character sum onto the appropriate even or mod 3 
representations in each case. 

The correlators in (|3.42 ) are also given by a simple contour integration. It can be shown 



that § 

1 d - ^{hfH-gys}^ 

\r\rr "\/ , *. t A \ \ * ± 



Ndt 



=- log X r hW x(A) = Y: -± fjr- (3.43) 



9 



1 f h (o) \ - — 7T\ 



where foj. is defined in (3Jj ). From the character expansion of section 2 we see that the 
left-hand side of (|3.43| ) arises from a derivative of the dual model free energy Sfl^[A,t*] with 



respect to the vertex weights t*, whereas the right-hand side is identical to (|3.7| ) with q — > —q. 
In the large- iV limit, we can therefore represent the dual model correlators in ( p.42|) by the 
contour integrations 

/ -(XA 3 f q ) ^-^-cf— e -^ H (h )+H h)) (344) 

Although it is possible to determine the complete set of observables of the dual matrix models 
in terms of the saddle-point solution at large- N, the situation is somewhat more complicated 
for the correlators of the Ai 3 and M.^ ensembles. The above discussion illustrates to what 
extent the large- iV saddle-point solution of the Itzykson-Di Francesco formula can be used to 
represent the observables of the matrix model. Those which do admit such a representation 
typically appear to be obtainable from derivatives of the free energy of the model and thus 
correspond in the random surface interpretation to insertions of marked loops on the surfaces. 
Thus, strictly speaking, the natural observable to compute in the large-iV limit of the 3- 
point matrix model is the free energy ( |3.19| ). In the next subsection we shall show how this 
calculation carries through. 



3.4.3 Free Energy 

The natural observable to compute in the triangulation model is the large- iV (genus zero) free 
energy (|3.19|) with I c = 1. Splitting up the integration range and setting pn(h) = 1 on [0,b], 
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we find 



S u = - [ dh p H {h)h [log(A 3 /i) - l] + - f a dh p H (h) [hlogh- (h - b) \og{h -b)-b\ 

3\ 6 2 - 3/2 



3 J b 



3 Jb 



1 r a r a b ( 3\ 6 — 3/2 

+ 6 1 Jb dk dk ' P H ^PH( h ') lo S \ h ~ ti\ + J { l °Z b ~ 2) + 4~ log A 

(3.45) 

The spectral density is found by computing the discontinuity ( |3.5| ) of the weight resolvent 
function ( 3.27 ) across the cut [b,a], which yields 

1 



Pn(h) 



71 



arctan 



2Jab(a -h)(h- b) 



(a + b)h - 2ab 



— 2 arctan 



(a-h)(h-b) 



h + yab 



, h G [6, a] 
(3.46) 

The double integral in the free energy ( |3.45| ) can be simplified by integrating up the 
saddle-point equation (|3.3|) for h G [b, a] to get 



-Tdh' PH (ti) 
Jb 



loxlh-til 



h 



1 - log(A 3 /i) + (h - b) log(h - b) - h log h + - log A 



+-f a dti PH (ti) \og(ti -b)+b \hg(X 3 b) 

J b L 



4 

b log b 



(3.47) 



Substituting ( |3.47| ) into ( |3.45| ) and integrating by parts, we find after some algebra 



S 



H 



1 r dh d PH {h) 



6 Jb 



dh 



h 



h 

h log h — — \og(h — b) + log(/i — 6) 

(/i> 



6(1 _-jlog(A-6) 



+-log6 1 + - log A (b + 15) + — log A 



6 



12 



6 



(3.48) 

where we have introduced the reduced weight average (h) = J b a dh pu{h)h = (h) — b 2 /2, and 
from ( |3.46|) we have 



dpn(h) h — 2\fab 



dh 7ry/(a-h)(h-b) 
The weight average (h) can be read off from (|3.27 ) to give 

(h) = -ri 2 /3 + r]-l/6 



h G [b, a] 



(3.49) 



(3.50) 



To evaluate the remaining logarithmic integrals in fl3.48| ), it is convenient to change vari- 
ables from h G [b,a] to x G [—1, 1] with h = |(a + b) + |(a — b)x. This leads to a series 
of integrals over elementary algebraic forms and forms involving logarithmic functions. The 
latter integrals can be computed by rewriting the integrations over x G [—1, 1] as contour inte- 
grals, blowing up the contours and then picking up the contributions from the discontinuities 
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across the cuts of the logarithms. The relevant integrals are then found to be 

log(l + rx) 



h(r) 
h{r) 

J(r) 



i 

dx — 

-1 y/T 



X' 



71 log 



1 + VT^r< 



1 x log(l + rx) 

dx 



i 



n -(i-VT 



1 , x 2 log(l + rx) 
dx 



i 



vT 



log(l + x) 



<ix . 

l (1 + rx) v 1 — x 2 



— 7T 



-1 



r 2 



dy 



(3.51) 



7T 



vT 



(i + n/)Vi -y 2 

f 1-r 

= 10K 



+ 



-log(g 



where < r < 1. Using the above identities and the boundary conditions ( |3.28| ), after some 
tedious algebra we arrive finally at the simple expression 



s„M = ilogA + ilog, f + ;L-g 



(3.52) 



for the free energy of the triangulation model, where we have ignored an overall numerical 
constant. 

To examine the scaling behaviour of the free energy about the critical point, we once again 
introduce a renormalized, continuum cosmological constant A and momentum II as in ( |3.37|) . 
From the boundary equation 3A 6 r/ 3 = r\ — 1 we can calculate derivatives of r\ with respect to 
the cosmological constant to get 

drj 3A^ 4 
dA T] — T) c 

which we note diverges at the critical point. From (|3.52| ) and ( |3.53| ) the first two derivatives 
of the free energy are then found to be 



(3.53) 



dS 



ii 



-Wfo-2)/6 



d 2 S 



(3.54) 



8A c ' v ' " ' OA 2 

Both derivatives in (|3.54 ) are finite at the critical point. Taking one more derivative yields a 
combination that does not cancel the singular (77 — ^ c ) _1 part of ( 3.53Q , i.e. there is a third 
order phase transition at the critical point. Substituting ( |3.37|) into the boundary equation 
(|3.28| ) gives II 2 ~ A near criticality A -> (see (|3.38| ),( ^39| )), and so substituting r] ~ r/ c -A 1/2 
into (|3.54|) and expanding about A = yields 



d 2 S 



H 



dA 2 



■ A 1 ' 2 + ... 



(3.55) 



which identifies the expected pure gravity string exponent (|3.41|) . 
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4 Incorporation of the Itzykson-Zuber Correlator 



In this section we shall give a more quantitative presentation of the evaluation of observables 
in the curvature matrix models and how they reproduce features of the triangulated random 
surface sum. This will also further demonstrate the validity of the weight splitting that was 
assumed to hold in the Itzykson-Di Francesco formula, i.e. that the leading order contributions 
to the partition function indeed do localize at N = oo onto the even representations of 
GL(N, C) that split symmetrically into mod 3 congruence classes. For this we shall present 
an explicit evaluation of the Itzykson-Zuber correlator J-'(h) using (|3.6| ). This will amount 
to an evaluation of the large-iV limit of the generalized Schur functions introduced in the 
previous section. 

Notice first that expanding both sides of ( |3.6|) in powers of q and equating the constant 
terms leads to the identity 

i = jf^ {n H {h) + m) (4-1) 

which, along with the normalization Jq dh ps{h) = 1 of the spectral density, implies that 

£dhJ r (h) = (4.2) 

Thus, at least for some range of the couplings £*, the function T{h) will be analytic in a 
neighbourhood of the cut of Tinih). This will be true so long as the Itzykson-Zuber integral 
does not undergo a phase transition which changes its analyticity features in the large- N 
limit. The equation ( |3.6|) can be used to determine T{h) once the coupling constants of the 
dynamical triangulation are specified. The simplest choice is a power law variation of the 
couplings as the coordination numbers are varied, 

qt* = t q - 1+p ■ t G R , p G Z (4.3) 

so that each vertex of the triangulated surface interacts with a strength proportional to t 
(with proportionality constant t p ) with each of its nearest neighbours. 

This simple choice of vertex weights in the triangulation model allows us to study more 
precisely how the curvature matrix model represents features of the random surface sum, and 
to further demonstrate that the actual saddle-point localization of the partition function is 
not onto some configuration of Young tableau weights other than the mod 3 representations. 
It will also allow us to examine how the observables of the model behave as the weights are 
varied in this simple case, without the complexities that would appear in the analytic form 
of the solution for other choices of coupling constants. For generic t ^ 1 (i.e. ^(h) ^ 0), the 
saddle-point solution of the matrix model must satisfy certain physical consistency conditions 
so that it really does represent the (continuum) genus zero contribution to the random surface 
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sum ( |1.3| ) with the choice of couplings ( |4.3| ). Using Euler's theorem V — E + F = 2 for a 
spherical topology and the triangulation relation 



J2 q = 2E = 3F (4.4) 

v*eG 3 

it follows that the planar surface sum in ( |1.3| ) is 

zP(X,t) = £ c G(0) (\H^) A[Gf)) (4.5) 



of " 



where the sum is over all planar fat-graphs of area A(G { 3 ] ) oc F(Gf ] ), and the constant 
c„(o) is independent of the coupling constants A and t. The perturbative expansion parameter 
is \ 6 t 2 ( p+1 ^ and the critical point, where a phase transition representing the continuum limit 
of the discretized random surface model takes place, is reached by tuning the expansion 



parameter to the radius of convergence of the power series ( [4.5|) . The critical line A c (i) thus 
obeys an equation of the form 

A c (t) 6 ■ t 2(p+1) = constant (4.6) 

The solution of the model for t ^ 1 should therefore have the same physical characteristics as 
that with t = 1, since in the random surface interpretation of the curvature matrix model the 
only effect of changing t is to rescale the cosmological constant of the random surface sum for 
t = 1 as A 3 — > A 3 t p+1 . This geometrical property should be reflected in the large- iV solution 
of the matrix model with a non-vanishing Itzykson-Zuber correlator. 

As discussed in ||, it is possible to invert the equation 

G(h) = e ^{h) + Hh) ( 4 _ 7 ) 

to obtain h as a function of G, by changing variables in the contour integral ( |3.6|) from h to 
G to get 

- L^g^ (48) 

where the contour Co encircles the cut [G(fe),0] with clockwise orientation in the complex 
G-plane. Then 

oo oo 

MG) = l + £^ + E^ (4.9) 
q =i G9 q =i 

where the coefficients 

determine the analytic part of the function h(G). The second equality in fl4.1U| ) follows from 
( 3.44 ) and the constant term in the Laurent series expansion Q4.9p is unity because of the 
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normalization of the spectral density pu ||- In the general case, the solution for G(h) as 
determined from ( |4.9| ) will be multi-valued. It was shown in || that the first set of sheets of 
the Riemann surface of this multi-valued function are connected along, and hence determined 
by, the cut structures of e^^, which map the point h = oo to G = 0. These sheets are found 
by inverting ( f4.9| ) with g q = 0. The remaining sheets are determined by the cuts of e nH ^ 
and are associated with the positive powers of G in (|4.9|). 

The solution is extremely simple, however, for the choice of couplings ( f4.3|) , as then the 
inversion of the equation ( f4.9| ) with g q = yields 

Gi(h)=t + 1 ^ (4.11) 

This solution has a simple pole of residue t p at h = 1, but no multivalued branch cut structure. 
Thus we expect that the singularities of will have a pole structure, rather than a cut 

structure. The remaining sheets of the function G(h) will be determined by the cut structure 
of e HH( - h \ They are attached to the "physical" sheet G\(h), on which the poles of and 
the cuts of e HH ^ lie, by these cuts. Note that the analytical structure of the character as 
determined by this part of the Laurent series in (^4.9|) is anticipated from the Schur character 
formula (|3~9p,(CT). 

With this observation we can in fact obtain a closed form expression for the Itzykson- 
Zuber correlator in the large- N limit, in contrast to the generic case || where in general one 
obtains only another discontinuity equation such as ( [373]) . The resolvent function (|3.4|) can 
be written as 

H H {h) = log [yzij + H H {h) (4.12) 

where 7tI#(/i) = J fe a dh! pH{h')/(h — h') is the reduced resolvent function associated with the 
non-trivial part of the density pn , and it has a branch cut on the interval [b, a] . Using ft4.2|) 
it can be written as 

Jc 2m h — h Jc G 2m h — h(G) 

where the contour C encircles the cut [b, a] of Tinih) and we have changed variables from 
h to G as above. The contour integral ( |4. in the complex G-plane can be evaluated for 
large-/i, which, by analytical continuation, determines it for all h. We can shrink the contour 
Cq down to an infinitesimal one Co hugging both sides of the cut [G(b), 0]. In doing so, from 
(4.S) it follows that we pick up contributions from the solution ( [4.11|) and the extra pole at 
G = t corresponding to the point h = oo. Thus integrating by parts we find that (|4.13| ) can 
be written as 

r dG id 1 1 

H H {h) = lo g G a (/0 ~^ t + f Co ^- \ qq [logGlog(fc - h{G))) - - log(/i - h(G)) | (4.14) 
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In the contour integral over Cq in ( fQjj ), the total derivative in G gives the discontinuity 
across the cut [G(b), 0], which is \og(h — b). The other term there evaluates the G = limit 
of \og(h — h(G)) determined by (JO) and ( |4.11| ). Using ( f4.12| ), we then have 



H H (h) = \ogG 1 (h) + \og 



h 



(h - l)t + tP 



(4.15) 



Since there is only the single (physical) sheet determined by the singularity structure of e^ h \ 
we have logGi(h) = ^(h) + 7i#(/i) on the physical sheet, and combined with ( [4.15 ) we arrive 
at the expression 

\h 



Hh) = log 



l)t + V 
h 



(4.16) 



for the large-iV limit of the Itzykson-Zuber correlator. As mentioned above, the fact that 
Q F{h) on iy a s j m ple pole of residue t p — t at h = is because there are no other sheets 
below Gi(h) connected by cuts of e^^. This is opposite to the situation that occurs in the 
Gaussian case (A = in (|1 . 1|) ) , where e HH ' fe ' has a simple pole of residue 1 at h — 1 and there 
are no upper branches above G\{h) connected by cuts of e Hff '^ [j|]. It can be easily verified, by 
blowing up the contour C, using the asymptotic boundary condition Tinih) ~ l/h+ 0(l/h 2 ) 
for large- /i and computing the residue at h = oo in ( p.6|) , that (|4.16|) consistently yields the 
weights ( |4.3|) . Furthermore, for p = t = 0, (|4.16|) reduces to the solution of H in the case 
where only the vertex weight t* is non-zero, while for t = 1 (i.e. A = A = 1), ( [4.1 6|) yields 
J-(h) = 0, as expected from its definition. 

Notice that for p / 1, F{h) here has a logarithmic branch cut between h = and 
h — i = 1 — t p-1 , so that strictly speaking the solution ( [4.1(j| ) is only valid for i < 0, i.e. 
0<t<lforp<0 and t > a for p > 1 (where its cut doesn't overlap with the cut [0, a]). 
Outside of this region the analytic structure of the Itzykson-Zuber correlator can be quite 
different. For p — 1, we have J-'(h) = logt and the only effect of the Itzykson-Zuber integral 
in the saddle-point equation ( |3.3|) is to rescale the cosmo logical constant as A 3 — > A 3 t 2 . This 
is expected directly from the original matrix integral ( |1.1|) , since for p = 1 the vertex weights 
can be represented as traces of the external matrix A — t • 1, while for p / 1 the GL(N, C) 
characters can only be defined via the generalized Schur functions fl3~9l ) , (|3.10|) . For p ^ 1, we 
shall now see that Q4.16 ) changes the analytic structure of the large- iV solution of the curvature 
matrix model, but that the t = 1 physical characteristics (the pure gravity continuum limit) 
are unchanged. 

The saddle-point equation ( |3 .3| ) with ( [4.16[ ) is then solved by replacing the t = 1 resolvent 
H„(h; A 3 , 1) in O by 



(h-a)(h-b) , f(s 

1 — log - — 

(s — a)(s — b) \ 



l)t + t p 



(4.17) 
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and compressing the closed contour C in ( 4.17 ) to the cut [t, 0] . Note that since i < the sign 
of the square root in ( [4.17| ) is negative along this cut. Working out the contour integration 
in ( }4.17| ) as before then gives 



U H (h; A 3 ,t) =H H (h;\H p+1 ,l) 

' h 2 (h(a + b)- 2ab - t(2h - a - b) - 2yJ(h - a)(h - b) (i - a)(i - b)) ' 



+ 2 log < 



(h - t) 2 (h(a + b)~ 2ab - 2y 'ab(h - a)(h -b)) 



(4.18) 

The endpoints a, b of the support of the spectral distribution function are found just as before 
and now the boundary conditions ( |3.28| ) are replaced by 



& = \H 2 ^rf t 



i-i=r- 1 = r h -zt t 



where 



« ~ 4 



-t-y/b 
and we have assumed that t ^ Q 



a-t+\b-t 



(4.19) 



(4.20) 



The boundary equations (|4.19|) are identical to those of the t — 1 model in ( |3.28| ) with 
the replacements ^ — > £ t = t 1 ^ p C, t , r\ — > fj t = t 1_p r]t and A 3 — > A 3 t p+1 . The weight average (/i) t 
corresponding to the l//i 2 coefficient of the asymptotic expansion of (|4.18|) is then 



(h)i 



+ t 



2(P-1) 




(4.21) 



where we have used (|4.19| ). Substituting ( |4.21| ) into (|3.33| ) yields (|3.35|) with r) —> fj t and an 
additional overall factor of t 2 ^" 1 ) which represents the overall factor in ( |4.5| ) (the linear in 
A term is from the Gaussian normalization of the partition function). Thus, although the 
precise analytical form of the solution is different, the critical behaviour (and also the Wick 
expansion) of the curvature matrix model for the choice of vertex weights ( |4.3| ) with t ^ 1 is 
the same as that for t = 1 . 

The correlators ( |3.32|) will all have a structure similar to those at t — 1, as in ( [1.21| ). The 

non-trivial analytical structure of the saddle-point solution for t ^ 1 is exemplified most in 

1 It can be easily seen that for p = t = the saddle-point equations are not satisfied anywhere. This simply 
reflects the fact that it is not possible to close a regular, planar triangular lattice on the sphere. As discussed 
in [|], U , in the matrix model formulations one needs to study "almost" flat planar diagrams in which positive 

curvature defects are introduced on the Feynman graphs to close that triangular lattice on the sphere. 
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the function ( |4.7| ), which in the present case is 

G{h;X 3 ,t) = G(h;X 3 t p+1 ,l) 

th 3 (h(a + b) - 2ab - t{2h - a - b) - 2sJ{h - a){h - b)(t- a){i- b)f 

X 2 

(h - i) 3 (h(a + b)- 2ab - 2^ab(h - a)(h - 6)) 

(h(a + b) - 2ab - t{2h - a - b) - 2yJ(h- a) (h -b)(t- a) (f - b)) 2 
A 3 tP(v^ + \/6) 2 (a -b){h- tf (h(a + b)- 2ab - 2y/ab{h - a)(h - bj) 

(4.22) 

The inverse function h(G) can be determined from ( [4.22 ) using the Lagrange inversion formula 



o \ fc-1 

a 



k=l 

where the function tp is defined from ( 4.22j ) by 



h(G) = G + j: B — ^(G) fc (4.23) 



^(/i) = /i - G(h) (4.24) 



Comparing with ( |4.9| ) for the choice of vertex weights ( |4.3| ), we can write down an expression 
for the generating function of the dual moments of the triangulation model (or equivalently 
for the 4-point and even-even models) 

oo /t > fP^n 00 1 / f) \ k ~ 1 



Because of the complicated structure of the function ( |4.22| ), it does not seem possible to 
write down a closed form for this generating function or systematic expressions for the dual 
moments. Nonetheless, ( |4.25| ) does represent a formal solution for the observables of the 



triangulation model. 

The above critical behaviour, when a non-trivial Itzykson-Zuber correlator is incorporated 
into the dynamical triangulation model, is anticipated from (|4.6| ) and thus yields a non-trivial 
verification of the assumptions that went into the derivation of the large-iV limit of the 
Itzykson-Di Francesco formula of section 2. The form of the function G(h;X 3 ,t) in fl4.22|) 



illustrates the analytical dependence of the saddle-point solution on the vertex weights t*. It 
also demonstrates how the analytical, non-perturbative properties of the random surface sum 
( |4.5| ) change at iV = oo, although the perturbative expansion of the free energy coincides 
with ( f4.5|) and the physical continuum limit ([4.6|) is unaltered. The discussion of this section 
of course also applies to the 4-point and even-even models with the appropriate redefinitions 
of coupling constants above, and also to the complex curvature matrix model where now the 
incorporation of the Itzykson-Zuber correlator using the saddle-point equation ( |3.22| ) leads 
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to the appropriate rescaling of the cosmological constant A 3//2 by £( p+1 )/ 4 as predicted from 
the graphical arguments of subsection 3.3 (see ( |3.24| )). The above derivation also suggests an 
approach to studying phase transitions in the large- N limit of the Itzykson-Zuber model, as 
it shows in this explicit example the region where the analytic structure of T{K) changes (i.e. 
t > 0) and consequently the region wherein a discontinuous change of the large- N solution 
appears. This could prove relevant to various other matrix models where the Itzykson-Zuber 
integral appears [0, [Llj] . The saddle-point solution above of the curvature matrix model can 
nonetheless be trivially analytically continued to all tel. This is expected since the random 
surface sum Q4.5|) is insensitive to a phase transition in the Itzykson-Zuber integral which 
appears in the large- iV solution of the matrix model only as a manifestation of the character 
sum representation of the discretized surface model. 



5 Complex Saddle Points of the Itzykson-Di Francesco 
Formula 

The curvature matrix models we have thus far studied have led to a unique, stable saddle- 
point solution at large- N. From the point of view of the Itzykson-Di Francesco formula of 
section 2, there is a crucial reason why this feature has occured, namely the cancellation 
of sign factors that appear in expressions such as ( |2.2| ). The models we have studied have 
been arranged so that there is an overall cancellation of such sign factors which arise from 
the splitting of the Young tableau weights into the appropriate congruence classes. When 
the weights are then assumed to distribute equally the resulting Vandermonde determinant 
factors stabilize the saddle-point and lead to a unique real-valued solution for the free energy 
and observables of the random matrix model. In this section we shall briefly discuss the 
problems which arise when trying to solve the matrix models when the sign variation terms 
in the character expansion formulas do not necessarily cancel each other out. 

The destabilization of the real saddle-point configuration of weights was pointed out in 
Q where it was shown that the configuration is complex for the Hermitian one-matrix model 
with Penner potential || 

Vp(XA) = -log(l-X) (5.1) 

in (|1 . 1| ) . The Itzykson-Di Francesco formula for this matrix model follows from replacing 
by X{h}{A) = X{7i}(l) oc A[h] in (|2.5| ), so that at iV = oo we have 

N/2 

Zf = c N \- N2 ' 4 Y, ^[h°] 2 A[h e ] 2 J] 0° - hf) e^*^)- 1 ] (5.2) 

h={h e ,h°} i,j=l 
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Now there is no problem with the decomposition of weights into the appropriate congruence 
classes, but, as we shall see below, the rapid sign changes of the mixed product over the even 
and odd weights destabilize the reality of the saddle-point configuration of Young tableau 
weights. In the previous models such mixed product factors did not pose any problem for the 
solution at large- iV because they appeared in the denominators of the character expansions 
and acted to make the more probable configurations of weights those with identical distribu- 
tions of even and odd weights, thus stabilizing the saddle-point at iV = oo. In ( |5.2|) , however, 
the mixed product Yli,j(h° — h^) appears in the numerator and thus acts to make the more 
probable configuration those with different distributions of even and odd weights. Thus when 
a symmetric distribution over h\ and h° in (|5.2| ) is assumed, this has the effect of destabilizing 



the saddle-point leading to a complex-valued solution. 

The matrix model with Penner potential ( [5.1| ) is equivalent to the standard Hermitian one- 
matrix model for pure gravity i.e. that with potential | tr X A in (|1 . 1| ) . Diagrammatically, 



a two-to-one correspondence between planar Feynman graphs of these two matrix models 
exists by splitting the 4-point vertices of the X 4 model as described in subsection 3.2 to give 
diagrams of the "even-log" model with potential — log(l — X 2 ) = — log(l — X) — log(l + X) 
(so that the face centers of the X A model are mapped onto the vertices and face centers of 
the even- log model). From the point of view of the Itzykson-Di Francesco formula, in the 
character expansion ( [2.5|) for the X A model we replace X s [h] by X^[h}. The resulting partition 
function Z 4 is a sum over mod 4 congruence classes of weights in which the distribution sums 
for the classes \i = 0, 2 and /i = 1, 3 decouple from each other (so that even and odd weights 
completely factorize) and each have the precise form at N = oo of the partition function (|5.2|) 
0, i.e. zf^ = (Zp^) 2 . This is just the graphical correspondence mentioned above. Thus 
the Itzykson-Di Francesco formula at least reproduces correct graphical equivalences in these 
cases. 

To see how the complex saddle-points arise in the character expansion of the Penner 



matrix model above, we assume that even and odd weights distribute symmetrically in (|5.2|) 
and define a distribution function pp(h) for the N/2 weights h\. Varying the effective action 
in (|5.2|) for the weights h e then leads to the large- iV saddle-point equation 

The corresponding resolvent function 7ip(h) can be determined by the contour integration in 
(|3.26| ) just as before using ( |5.3j ) and we find after some algebra 

^hLoJ'^-^- 6 ' + ->(*-»»' } (5 . 4) 

3 [ A (h(a + b)- 2ab + 2y/ab(h - a)(h - b)] 
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Expanding ( |5.4j ) for large- h then leads to the boundary conditions 

£=jjfa-l) , @ 3 AV = (^-1) 3 (5.5) 

Consider the structure of the solutions to the boundary conditions G5.5|) . Again, the Wick 
expansion of the original matrix integral is analytic about A = 0, and it can be shown that 
the free energy is analytic about r](X = 0) = 1. We should therefore choose the branches of 
the quintic equation in ( |5.5| ) which are regular at A = 0. There are three solutions which obey 
this analyticity condition and they are given by the iterative relations 

Vn = l + l^X 2/3 (Vn? /3 , n = 0,1,2 (5.6) 

where G Z3 is a non-trivial cube root of unity. The remaining two solutions behave near 
A = as 77 ~ ±A _1 . The discrete Z^-symmetry of the regular saddle-point solutions (5.6) 



seems to be related to the fact that the Schwinger-Dyson field equations of this matrix model 
determine the function Gp(h) = e Hp ^ as the solution of a third-order algebraic equation Q 



\h 3 Gp - \h 2 (l + h)G 2 P +\\-h + (l - Vl - 12A) (1 - 12A) 

y 64oA v 



hG P + h 2 = (5.7) 



Initially, the endpoints a, b of the support of the spectral density lie on the positive real 
axis, so that one expects that only the real branch 770 is a valid solution of the matrix model. 
However, the perturbative expansion parameter of the free energy Sp(r)o) would then be A 2 / 3 . 
It is easy to see by a Wick expansion of the original matrix integral that the genus zero 
expansion parameter is in fact A 2 . Furthermore, one can analyse the analytic properties of 
the solutions to the quintic boundary equation in (|5.5|), and even determine the critical point 
A c which in this case is the point where the two real and positive roots coalesce. For A > A c all 
three roots which are analytic about A = become complex-valued. This critical behaviour 
is similar to that discussed for the cubic boundary equation (|3.28| ) in subsection 3.4, and so 



apparently lies in the same universality class as the earlier models. However, the critical value 
A c determined this way does not agree with known results [1], so that the free energy Sp 
determined from the character expansion does not count the Feynman graphs of the Penner 
model correctly. 

The structure of complex saddle-points arises in many other matrix models. For example, 
for the quartic-type Penner potential 

vj?\XA) = -log(l-X 4 ) (5.8) 

the boundary conditions determining the endpoints a, b of the support of the distribution 
function are 

£ = f (V - 1) , (£) AV = {V~ I)' (5-9) 
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Again, there are three regular solutions of ( |5.9| ) at A = and the real root leads to an 
expansion in A 4 / 3 , whereas the Wick expansion can be explicitly carried out and one finds 
that the perturbative expansion parameter is A 4 . 

All the models we have studied for which the matrix model saddle point does not reproduce 
the graphical expansion share the feature that the constraint of regularity at A = yields 
multiple solutions of the endpoint equations. Choosing the real root (or any other single root) 
leads to the wrong solution of the matrix model. Thus it appears that the saddle-point should 
be defined by some kind of analytical continuation that extends the support of the spectral 
density pp into the complex plane and takes proper account of the multiple root structure. It 
would be interesting to resolve these problems and determine the general technique for dealing 
with such complex saddle-points. It would also be interesting to discover any connection 
between this saddle-point destabilization and the well-known occurence of (single-branch) 
complex saddle points for the eigenvalue distributions in generalized Penner models ||16|| . 



6 Conclusions 

In this paper we have shown that the character expansion techniques developed in [[|, [5] can be 
applied to odd potentials. We have demonstrated that the splitting of weights into congruence 
classes other than those of the even representations of GL(N,C) leads to a proper solution 
of the matrix model, provided that one writes the weight distribution that appears in the 
character expansion over the appropriate congruence elements. The Itzykson-Di Francesco 
formula then correctly reproduces relations between different models. 

From a mathematical perspective, the results of this paper raise some questions concern- 
ing the large- N limit of the Itzykson-Di Francesco formula. For instance, a random surface 
model with a discrete Z p -symmetry corresponding to, say, a p-colouring of its graphs will be 
described by a curvature matrix model with the Z p -symmetry A — > uj p A, uo p e Z p . This sym- 
metry will be reflected in the Itzykson-Di Francesco expansion as a localization of the group 
character sum onto mod p congruence class representations of GL(N,C). The appropriate 
solution of the model at N = oo will then involve resumming the Young tableau weights 
over the mod p congruence classes and assuming that the even-odd decomposition factorizes 
symmetrically over these classes. However, it is not immediately clear why such a symmetry 
assumption of the character expansion at large- N gives the appropriate solution of the dis- 
cretized planar surface theory (although a mapping onto a complex matrix model indicates 
how this should work). At this stage there seems to be a mysterious "hidden" symmetry at 
play which makes the large-iV group theoretical approach to solving these random surface 
models work. Furthermore, other intriguing features of the Itzykson-Di Francesco formula, 
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Figure 5: An example of a regular fractal-type graph that appears in the dynamical triangulation. 



such as the appearence of phase transitions in the Itzykson-Zuber integral represented through 
the large- iV limit of generalized Schur functions, are purely a result of the character expan- 
sion representation and correctly conspire to yield the proper solution of the random surface 
model. It would be interesting to put all of these features into some systematic framework 
for dealing with curvature matrix models in general. 

From a physical point of view, there are many random surface models which are best dealt 
with using the dynamical triangulations studied in this paper, and it would be interesting to 
exploit the relationships with the even coordination number models to study the properties of 
these theories. For instance, one sub-ensemble of the random surface sum ( p..3[ ) is the collection 
of regular fractal-like graphs (Fig. 5) which were shown in [|17[] to dominate, in the high- 
temperature limit, the surface sum for two-dimensional quantum gravity coupled to a large 
number of Ising spins when restricted to two-particle irreducible Feynman diagrams. These 
fractal graphs can be characterized by 3-point graphs G3 where only the dual coordination 
numbers 

q* k = 3-2 k , k>0 (6.1) 

are permitted. Here k is the order of the fractal construction obtained inductively by replacing 
each 3-point vertex of an order k — 1 fractal graph with a triangle^. The ensemble of fractal 
graphs corresponds to a branched polymer phase of two-dimensional quantum gravity || . We 
have shown that when the 3m-sided polygons in ( |1.3|) are weighted with a power law variation 
with the number of sides, the continuum limit of the model lies in the pure gravitational phase. 



2 Note that the total number of triangle sides along each outer side of the fractal-like structure of a single, 
order k + 1 fractal graph is 2 k — 1. Thus to be able to close the set of 3-point graphs of dual coordination 



numbers (p.l| ) on a spherical topology (corresponding to N — 00 in the matrix model), one needs to glue an 
order k and an order k + 1 fractal graph together along their three external corner legs (see Fig. 5) so that 



the valence of the dual vertices of the faces joining the two fractal structures will coincide with (6.1 
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The curvature matrix model ( |1 . 1| ) with the dual vertex weights arranged as discussed in this 
paper can thus serve as an explicit model for the transition from a theory of pure random 
surfaces (associated with central charge D < 1) to a model of branched polymers (associated 
with D > 1). This might help in locating the critical dimension D c > 1 where the precise 
branched polymer transition takes place. 
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